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ABSTRACT 

T Ceti (HD10700), a G8 dwarf with mass 0.78 M©, is a close (3.65 pc) sun-like 
star where 5 possibly terrestrial planet candidates (minimum masses of 2, 3.1, 3.5, 
4.3, and 6.7 M®) have recently been discovered. We report abundances of 23 elements 
using spectra from the MIKE spectrograph on Magellan. We find [Fe/H] = -0.49 
and Teff = 5387 K. Using stellar models with the abundances determined here, we 
calculate the position of the classical habitable zone with time. At the current best fit 
age, 7.63;|;°-5’ Gy, up to two planets (e and f) may be in the habitable zone, depending 
on atmospheric properties. The Mg/Si ratio of the star is found to be 2.01, which is 
much greater than for Earth (~1.2). With a system that has such an excess of Mg to 
Si ratio it is possible that the mineralogical make-up of planets around r Ceti could 
be significantly different from that of Earth, with possible oversaturation of MgO, 
resulting in an increase in the content of olivine and ferropericlase compared with 
Earth. The increase in MgO would have a drastic impact on the rheology of the mantles 
of the planets around r Ceti. 

Subject headings: astrobiology, planets and satellites: composition, planets and satel¬ 
lites: interiors, stars: abundances, stars: individual(T Ceti.) 


1. Introduction 


With the number of extrasolar planets increasing rapidly, it appears that Earth-sized planets 
in their host star’s habitable zone (HZ, (e.g. Kasting, Whitmire, & Reynolds||1993[ [Kopparapu et 
al.|2013a| )) should be numerous ( |Marcy et al.|2000[ Gaidos|2()T^ pressing & Charbonneau|20T3 


Kasting et al.||2013| [Petigura et al.||20l3| ). In two decades we have progressed from having no 
candidate planets to having too many to practically search for detectable biosignatures. A more 
nuanced analysis than the location of a planet relative to the instantaneous HZ is necessary to 
choose among candidates. One way of conceptualizing this process is through a “detectability 
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index” (DI) would account for the ability of a planet to host life, its ability to maintain biosigna¬ 
tures, and the observational teehniques that determine what qualifies as a detectable level. The first 
two terms depend in part upon the characteristies of the host star that determine the position and 
time evolution of the HZ. Earth’s biosphere took ~ 2 Gy to produce potentially detectable ehanges 
in the amount of non-equilibrium species in the atmosphere. Hence a planet that entered the HZ 
only a few hundred million years ago as the star gradually increased in luminosity may well be 
habitable and even inhabited, but not have deteetable biosignatures. Both the stellar and planetary 
properties that determine habitability and detectability for terrestrial planets depend fundamentally 
on ehemieal eomposition. 


The determination of elemental abundances has led to the knowledge that every star has its 
own unique fingerprint. Variations on element to element seales exist between all stars for every 
element. As a community, we have begun the exploration of extreme cases, especially for the 
significant rock- forming elements Mg, Si, O and C (Delgado Mena et al.|2010[ Bond et al.|2010b 


Unterbom et al.|20T3 ). Sinee the stellar enrichments are likely primordial, it should lead to planets 
forming from the same enriched material Santos et al. (2001[ [2003); Valenti & Fischer (2005); 


Bond et al. (2006). 


Although stellar abundanees do not eorrelate direetly with planetary abundances, planet as¬ 


sembly simulations have been used to estimate the final planetary compositions (OBrien et al 
2006t[Bond et al.|2010b| ). [Bond et d^ ( |2010b] ) found that substantial stellar enriehments eould lead 
to planets with compositions enhanced relative to Fe in partieular elements, espeeially for the 4 
major roek-forming elements Si, Mg, C, and O. This diversity in bulk eompositions could lead to 


variation in mantles and core eompositions of terrestrial planets (Elkins-Tanton & Seager 2008). 


Reeently mueh attention has been given to earbon-rich systems (Bond et al. 2010b Kuchner & 
Seager|2005| [Madhusudhan et al.||2012| ), with [Bond et al.| ( [2010b[ ) finding 21 of their 60 stars hav¬ 
ing C/0 ratio above 0.80, believed to be a threshold for earbon and earbide dominated planets, 
as opposed to solar system-like silieate dominated planets. There are few observations of stars 


with extreme Mg/Si, with only a handful of stars over Mg/Si>1.9. Bond et al. (2010b) used initial 
values for HD177830 (Mg/Si of 1.93) in simulations, finding an olivine-rieh planet outside of 0.3 
AU, with up to 22% wt Mg and Mg/Si = 1.71, mueh greater than that of Earth. 


T Ceti is of long-standing seientific and popular interest as a nearby (3.65 pc) G8 dwarf. The 
recent discovery of planet candidates around r Ceti ( |Tuomi et al. ||2013| ) gives us an example of a 
nearby sun-like star with Super-Earth terrestrial planets. Two of these planets, r Ceti e and f, eould 


be in the HZ. These planets are estimated to have masses of 4.29 + 2.00 and 6.67 + 3.50M® (Tuomi 


et al. 1 2013). A high profile and potentially promising system sueh as r Ceti warrants a careful 


analysis. In this study we evaluate properties of the system, as a simple eoneeptual example and 
how it’s properties eould affeet habitability and the DI from a ehemieal perspeetive, eontributing 














































-3- 


to a complete and accurate characterization of the system. 

Section|^describes new observations from the MIKE spectrograph on Magellan used to derive 
abundances for 23 elements. Section [^discusses the possible effects of Mg-rich compositions on 
the internal dynamics and structure of terrestrial planets. New stellar models and HZ predictions 
using the newly determined abundances are presented in Section]^ The possible implications for 
the habitability of r Ceti e and f are discussed in Section]^ 


2. Abundance Determinations 


2.1. Observations 


The spectrum for t Ceti was taken by Paul Butler from the Carnegie Institute of Washington 


as part of the Magellan Planet Search Program, using the MIKE echelle spectrograph!Bernstein et 
^|2003 ) on the 6.5 m Magellan II telescope. The spectrum has R 50,000 in the range 4700 A to 
7100 A. Only the red chip data was used for abundance determination. Wavelength calibration was 
carried out by iodine absorption cell( |Marcy & Butl^|1992| ). The range of S/N for the sample is 
70 to 800. We took the uncalibrated spectra through the abundance finding process, starting with 
wavelength calibration in IDE and combining all of the orders into a single EITS file. Equivalent 
widths were measured by hand in IRAK Stellar parameters were determined using an iterative 
routine that simultaneously optimizes stellar effective temperature T'eff, surface gravity log(g), and 
micro turbulent velocity, using balance constraints of the abundances derived from neutral (Ee I) 
and ionized (Ee II) iron. The script started at solar values (effective temperature T^ff = 5777 K, 
surface gravity log(g) = 4.44, microturbulent velocity ^ = 1 km/s, and [Ee/H] = 0) in the ATEAS9 
model atmosphere grid. The stellar '^eff and ^ were determined through balancing the reduced 
equivalent widths and excitation potentials for Ee I lines. This process yielded an [Ee/H] for the Ee 
I lines. The [Ee/H] from Ee II is primarily influenced by the stellar surface gravity and was forced 
to match the [Ee/H] from Ee I lines by adjusting the surface gravity in 0.01 dex intervals until 
the mean abundance of Ee I lines matches that of Ee II lines. Uncertainties in temperature were 
determined by adjusting the temperature solution until the correlation between the abundances and 
excitation potential reached a Icr linear correlation coefficient for the given number of lines. This 
was done by slowly increasing the temperature until there was a linear Icr deviation from the calcu¬ 
lated value. The same Icr approach was used for the ^ uncertainty but with the abundances versus 
reduced equivalent widths. The errors for [Ee/H] were determined by using a standard deviation 
of all the line abundances for Ee I, combined with the change resulting from varying the Tspec and 
^ by their Icr uncertainties. The elemental abundances were determined using the spectral anal¬ 


ysis program MOOG (Sneden 19731, with a model atmosphere created by the program mspawn 
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through ATLAS9 ( Kurucz||1993 1. The errors for these elements are determined from the standard 
deviation of abundances given for each spectral line. At such low metallicities, the measurement 
of oxygen for the weak lines found in this spectrum becomes difficult, hut remain impor¬ 


tant. Therefore, the errors for oxygen were obtained using methods from|Bertan de Lis et al. 


( |2015| ). These are a combination of the differences in abundances between the spectral lines, 
by varying the stellar parameters by the parameter errors( 53 K, .06 dex surface gravity, 
0.08 dex metallicity, and 0.23 kms~ 1 microturbulent velocity) to account for the sensitivity 
to these parameters, and by uncertainties in the equivalent widths due to continuum place¬ 
ment. Since the oxygen lines at this low of a metallicity have such shallow spectral features, 
the placement of the continuum yields the largest source of error. 


2.2. Abundances 


We find T^ff = 5373K + 53 K and log (g) = 4.55. r Ceti is relatively metal-poor, with 
[Fe/H] = -0.49 + 0.08. Table gives stellar parameters and abundances [X/H]. Values for the 
absolute number ratios C/0 and Mg/Si are also provided. These ratios are taken before normalizing 
to any solar spectrum. Due to the low metallicity and the spectral range, the C/0 ratio has a fairly 
large uncertainty. Table gives the line list and line parameters. The solar values were calculated 
using a solar spectrum and compared to Tau Ceti on a line to line basis to calculate the solar 
normalized values in but no solar values were used when dealing with elemental ratios such as 
Mg/Si. 


This work yields Mg/Si = 1.78 + 0.04, much higher than Earth Mg/Si =1.21 ( [Kargel & Lewis 
1993] ) and the sun’s Mg/Si = 1.05 ( |Asplund et al.||2005[ ). The absolute abundance of Mg is higher 
than the sun, despite the star’s low metallicity. Ca is also supersolar. The Al/Fe is enhanced, but 
the absolute abundance of Al is slightly less than solar. 

The values of Mg/Si in the literature vary considerably. Tau Ceti appears to have a high Mg/Si 


ratio compared to the sun (1.78 in this study, with a mean of 1.48 in literature sources!Thevenin 


1998[|Bond et al.|2006[|2008||Takeda|2007[ [Neves et al.|2009^|Adibekyan et al.|2012[ ). The values 


range from 1.29(Adibekyan et al.|2012) to 1.78. However, since the Mg/Si ratio is dependent on 2 


separate unnormalized values, small deviations in both values can lead to larger differences in the 
ratio. It can be more useful to check the Mg and Si values individually to see the true variation 
in the Mg/Si ratio. The absolute (not solar normalized) values for Si from 9 different studies. 


including this one 

(Thevenin|1998[ Bodaghee et al.|2003t Valenti & Fischer|2005 

;Gilli etal.|2006; 

Bond et al.|2006; 

Takeda|2007 

Neves et al.|2009| Adibekyan et al.|2012 

) range from 7.12(Gilli et 

al.|2006|to7.33( 

Theveninll998), with a mean of 7.22 and median of 7.2. The value reported here 


is exactly the median of 7.2. The Mg values from 6 different studies, including this one(Thevenin 
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1998t |Takeda|[2007t [Neves et aL]|2009[ [Milone et al.||2011[ [Adibekyan et al.||2012| ) range from 7.3 


to 7.53 with a mean of 7.41 and a median of 7.43. This study finds a Mg/H value of 7.45, whieh 
onee again is elose to the median value. Varying the our values by the errors from this study (+0.05 
for Mg and ±0.04 for Si) the value for Mg/Si ean be decreased to as low as 1.33 or as high as 2.19. 
Both this work and the median values of other careful determinations in the literature point to a 
high Mg/Si ratio for r Ceti, and even the lowest reported values are still substantially higher than 
the sun. Even if r Ceti eventually proves to be at the low end of this range, its example indicates 
that the possibility of exoplanets with exotic silicate mineralogy needs to be explored. 


3. Internal Structure and Dynamics of Mg-rich Terrestrial Planets 


Though the range of Mg and Si from different groups vary more than quoted errors for this 
star and many others! Hinkel et al.|2014 ), r Cetis Mg/Si is significantly higher than the solar value 
of 1.05 (Asplund 2005). Regardless of the precise value, the example of r Ceti encourages us 
to explore the geophysical effects of distinctly non-terrestrial compositions on extrasolar rocky 
worlds. The theoretical behaviors of high Mg/Si planets are not based on specific cutoffs, but 
rather a gradual change from Earth-like to Mg-rich. 

The estimated sizes of the candidate exoplanets span a considerable range, from 2.0 to 6.6 M®. 
It is considered likely that they qualify as terrestrial planets, but it is uncertain without an estimate 
of bulk densities. The following analysis assumes that the planets are terrestrial Super-Earths or 
that smaller, undiscovered terrestrial planets are present. It also assumes that Mg/Si in terrestrial 
planets scales roughly with Mg/Si in the star. Stellar abundances most likely represent primordial 
compositionf jSantos et aH200T| [2003 1 [2005^ [Valenti & Eischer||2005[ [Bond et aH2006[ ) . We know 
from the example of Earth that there is not a 1:1 correspondence between the abundance ratios of 
star and planet (see Section 2.2| ), but for our solar system at least, the numbers are similar, with 
the terrestrial value being higher than the solar. This work speculates on the effects of an extreme 
Mg/Si ratio such as that found for r Ceti by some groups, including this study, on terrestrial planets. 
This kind of analysis would benefit from tighter constraints on the planetary compositions, but 
stellar properties can provide useful boundary conditions. Such constraints would likely come 
from models of planet assembly from a protoplanetary disk of r Ceti’s primordial composition. 
The compositions of planets can differ from a star’s due to various processes occurring during 
planet formation and evolution. Eor example, the Mg/Si ratio of Earth’s upper mantle is greater by 


20% than that of our Sun’s photosphere. Bond et al. (2010a I indicates that stellar enhancements 


will be reflected in planets, even if the exact ratio is not the same. While solubility of Mg may be 
very small for the metallic iron core and some amount of Si can exist as an impurity in the core, 
most Mg and Si may remain in the rocky mantle. 
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Table 1: Abundances of r Ceti 



Value 

± 

Temp 

5373 

53 

MtVel 

1.30 

0.23 

log(g) 

4.55 

0.06 

[Fe/H] 

-0.49 

0.08 

[C/H] 

-0.84 

0.14 

[0/H] 

-0.48 

0.12 

[Na/H] 

-0.24 

0.06 

[Mg/H] 

0.14 

0.03 

[Al/H] 

-0.07 

0.01 

[Si/H] 

-0.28 

0.04 

[S/H] 

-0.49 

0.11 

[Ca/H] 

0.04 

0.10 

[Sc/H] 

-0.32 

0.09 

[Ti/H] 

0.16 

0.08 

[V/H] 

0.16 

0.05 

[Cr/H] 

-0.33 

0.22 

[Mn/H] 

-0.26 

0.18 

[Co/H] 

-0.21 

0.04 

[Ni/H] 

-0.21 

0.12 

[Cu/H] 

-0.20 

0.15 

[Zn/H] 

-0.83 


[Y/H] 

-0.26 

0.30 

[Mo/H] 

-0.18 


[Ba/H] 

-0.44 

0.02 

[La/H] 

-0.16 


[Ce/H] 

0.23 


[Nd/H] 

-0.03 

0.38 

[Eu/H] 

-0.11 


C/0 

0.27 

0.21 

Mg/Si 

1.78 

0.06 
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The mantles of planets with Mg/Si approaehing 2 will be almost entirely made of olivine, 
(Mg,Fe) 2 Si 04 , with very little or no pyroxene, (Mg,Fe)(Al,Si) 03 . This is in sharp eontrast with the 
mineralogy of Earth’s mantle where the proportions of these two minerals are similar to eaeh other. 
Furthermore, the Mg-rieh mantle eould even contain ferropericlase, (Mg,Fe)0, at shallow depths. 
Ferropericlase exists only in the lower mantle (700-3000 km depths) in Earth. If the Mg-rich 
terrestrial exoplanets have surface volcanism through (partial) melting of the deep mantle similar 
to Earth, such as some hot spots and mid oceanic ridges, the high Mg/Si ratio will result in very 
different magnitude of melting and composition of melts. Mantle melting processes are believed 
to be responsible for the formation of the crust on Earth. 


If the composition remains similar throughout the mantle, the deeper mantle may be com¬ 
posed mainly of Mg-silicate perovskite (Mg,Ee)(Al,Si )03 and ferropericlase in similar proportion, 
unlike the Earth’s lower mantle where Mg-silicate perovskite is dominant (70%). It has been pre¬ 
dicted that Mg-silicate perovskite may break down into oxides (ferropericlase and silica) at 1 TPa 
(Umemoto et al. 20061. However, the known terrestrial planet candidates around r Ceti may not 
have sufficient mass to reach this pressure. Therefore, except for the shallow depths (< 600 km), 
the mantles of these large terrestrial planet candidates should be composed mainly of the phases of 
the Earth’s lower mantle but with different proportions. 


The increased amount of ferropericlase may have a profound impact on mantle dynamics. Eor 
the lower mantle, an increase in ferropericlase expected from the elevated Mg/Si ratio (- 2) could 
decrease the viscosity of the lower mantle by two orders of magnitude ( Ammann et al.|20TT ). The 
severe decrease in mantle viscosity could result in more vigorous convection, potentially affecting 
surface tectonics. Even if improved determinations of Mg/Si in r Ceti converge on a lower value 
than that found in this study, such novel geophysics are likely to appear in other exoplanet systems. 


4. Habitable Zone 


We evaluate the evolution of HZ boundaries defined by theoretical stellar evolution tracks 
coupled with the radiative-convective atmospheric models of [Kopparapu et ah ( 2013a|b ). Evolu¬ 
tionary tracks were created with the TYCHO code ( Young & Amett|2005 ). TYCHO is a ID stellar 
evolution code with a hydrodynamic formulation of the stellar evolution equations. It uses OPAE 
opacities ( |Iglesias & Rogers||I996^ [Alexander & Eerguson||I994] [Rogers & Nayfonov||2002[ ), a 
combined OPAE and Timmes equation of state (HEEMHOETZ) ( Timmes & Arnett] 1999 [[Rogers 


& Nayfonov 20021, gravitational settling (diffusion) (Thoul, Bahcall, & Eoeb 19941, automatic 
rezoning, and an adaptable nuclear reaction network with a sparse solver. A 177 element network 
terminating at ^"^Ge is used throughout the evolution. The network uses the latest REACEIB rates 


( Rauscher & Thielemann |2()00| Angulo et al.|1999 Iliadis et al.|2001 Wiescher et al.|20()6 ), weak 
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rates from Langanke & Martmez-Pinedo (20001, and screening from Graboske et al. (19731. Neu 


trino eooling from plasma proeesses and the Urea proeess is ineluded. Mass loss is ineluded but is 


trivial for a IM© main sequenee star. (Heightened early mass loss seen in some young stars (Wood 
et al.|2005) is not included. This may play a role in early HZ evolution and be relevant to the “faint 


young sun” paradox). It incorporates a deseription of turbulent convection ( 

Meakin & Arnett 2007 

Arnett, Meakin, & Young|2009 

Arnett, Meakin & Young||2010^1 Arnett & Meakin||2011 1 which is 


based on three dimensional, well-resolved simulations of convection sandwiched between stable 
layers, whieh were analyzed in detail using a Reynolds deeomposition into average and fluctuating 
quantities. It has no free parameters to adjust, unlike mixing-length theory (MLT). 


We used an initial mass of 0.784 M© (Tang & Gai 2011). The initial eomposition has the 


elemental ratios derived in Seetion 2.2, sealed so that when diffusion is ineluded the modeled 


surfaee composition matches the measured abundances at the best-fit model age. Traeks were 
computed to the end of its main sequenee lifetime. In order to find the best-fit age, we step through 
the evolutionary traek and evaluate: 


2 i.^mod^obs) (.^mod^obs} 

X = -^^- 


cr 


R 


(Tr 


( 1 ) 


where Rmod and Robs are the model and observed radii, Lmod and Lobs are the model and observed 
luminosities, and cr^ and ctl are the observational errors in radius and luminosity. The uneertainty 
in age is given by the ages at whieh the evolutionary traek erosses the observational error box. The 
adopted radius is 0.78 HR© , whieh is the mean of the interferometric and asteroseismological radii 
(0.773+0.024Ro and 0.07916+0.016Ro , respectively) determined by ( Tang & Gai|2011 ). For our 
error in R we use the error of the interferometrically defined radius, whieh eneompasses the range 
of uncertainty of the asteroseismologieal value. This provides us with a conservative estimate of 
the age uncertainty. We adopt L = 0.504Lo ( [Teixiera, et al.|[2009| |Pijpers|[2003[ [Pijpers et al. 
2003| ). [Tang & Gai] ( |2011| ) find log g - 4.53 using their asteroseismologieal analysis, similar to our 


speetroseopioally-determined log g = 4.55. The age of the system from this analysis is 7.63;|;j ^ Gy. 

We evaluate the elassieal HZ around r Ceti, where liquid water is stable at the planetary sur¬ 
faee given reasonable assumptions about the planetary atmosphere. Using the TYCHO simulation, 
we ean prediet the loeation of the inner and outer edges of the HZ as a funetion of time, and ask the 
question of whether a planet would be in the HZ long enough to harbor potential life. Locations 
for the inner and outer edges of the HZ are ealeulated using the stellar luminosity and temperature 
parameterizations from |Kopparapu et'aL](|2013a|bl). [Kopparapu et al.|(|2013a|) provide updated 1- 


D radiative-convective elimate models (originally from Kasting, Whitmire, & Reynolds (1993)) to 


obtain new estimates for HZ widths around F, G, K, and M-type stars. Rayleigh seattering by water 
vapor in the atmosphere of exoplanets with new seattering eoeffieients is ineluded. The models as¬ 
sume Earth-mass planets with water- (eorresponding to the inner HZ edge) or CO 2 (corresponding 
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to the outer HZ edge) dominated atmosphere. Both inner and outer edge ealeulations are based on 
’’inverse elimate modeling” where the planetary surfaee temperature is speeified, and the models 
are used to determine the stellar flux needed to sustain that given temperature. One potential issue 
is that the r Ceti planet candidates are all super-Earths. Models of planets with higher masses than 
earth and thicker atmospheres are discussed, but parameterizations are not provided. For the cases 
discussed, the edge of the HZ moves by at most a few percent in radius, so for the present we adopt 
the|Kopparapu et aL](|2013a|) parameterizations. 


Figure shows the inner and outer radii of the HZ for the cases in Table 3 of 
et al.| ( |2013a|b[ ). The orbits of the putative planets are indicated by the horizontal dashed lines. 
The vertical line is the best-fit age of the system at 7.63 Gy with the age uncertainty indicated 
by the gray box. The solid lines show the inner edge of the HZ for, in order of distance from 
the star, the Recent Venus, Runaway Greenhouse, and Moist Greenhouse cases. The outer edge 
is shown by the dotted lines for the Maximum Greenhouse and Early Mars cases. On the main 
sequence the radii of the HZ limits increase monotonically with time as the stellar luminosity 
increases. The three inner candidates are interior to the HZ even for the most optimistic prediction. 
Candidates e and f are potentially interesting marginal cases. For the Recent Venus case (assuming 
that Venus moved out of its HZ 1 Gy ago) planet e is reaching the end of its habitable lifetime. At 
the best-fit age of r Ceti it is right at the edge of the HZ. Future direct observations from space- 
based interferometers/occulters of this world have the potential to place constraints on climate 
models. For the theoretically-motivated greenhouse cases, candidate e is always interior to the 
HZ. Candidate f is potentially more promising, at least judging by the theoretical climate models. 
The outer edge of the HZ is approaching the position of planet f. Given the uncertainty in the 
stellar age, arising largely from the observational uncertainty in the star’s luminosity, candidate 
f is between 0 and 1.5 Gy from entering the HZ defined by the Early Mars case (assuming Mars 
was habitable 3.8 Gy ago). The range is extended by ~1 Gy for the Maximum Greenhouse case. 
Even in the most pessimistic case, the planet will have about 7 Gy of habitable lifetime until the 
end of the main sequence, plus additional time while the star traverses the subgiant branch. From 
a detectability standpoint, however, f is a poor candidate. At best, the planet has been in the HZ 
for <<lGy under these assumptions. The rate of change of L and T^ff as a function of time means 
that cases similar to f where a planet enters the HZ in the latter part of the star’s life are more 
common than planets that have been in the HZ since early times. This serves as a reminder that 
the present “habitability” of a planet does not necessarily indicate that it is a good candidate for 
detecting biosignatures. The temporal evolution of the system must be taken into account. 


Kopparapu 


The general characteristics of the r Ceti HZ are also interesting in the event that other planets, 
perhaps closer to the mass of Earth, are discovered. Table summarizes the range of orbits with 
given durations of habitability, assuming a stellar age of 7.63 Gy. We find ranges for orbits that are 
continuously habitable for the entire main sequence, orbits that remain habitable for at least 4 Gy 
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(the length of time required for metazoan life to develop on earth), and orbits that are eurrently 
habitable and have been so for at least 2 Gy. This last ease refleets the approximate amount of time 
taken for life to make potentially deteetable ehanges in the eomposition of Earth’s atmosphere (in 
partieular the presenee of free oxygen, methane, and nitrous oxide). Two eases are tabulated, a eon- 
servative ease for the minimum size HZ using the Moist Greenhouse and Maximum Greenhouse 
preseriptions, and an optimistie ease using the Reeent Venus and Early Mars elimate models. Eor 
the eonservative ease there are no eontinuously habitable orbits. In this dire situation the maximum 
duration of habitability is a mere 16 Gy, sinee a 0.78© star has a very long lifetime. Of the known 
eandidates, the best ease for a planet whieh has deteetable atmospherie signatures of an aetive 
biota is r Ceti e. The most favorable loeation for a hypothetieal sixth planet that has developed 
atmospherie biosignatures is the 0.45 AU range between 0.73 and 1.18 AU. 


5. Conclusions 

The r Ceti system has been the objeet of speeulation about extraterrestrial life in both the 
popular imagination and the seientifie literature for deeades due to its proximity to Earth and the 
star’s sun-like eharaeteristies. The reeent diseovery of Super-Earth mass planet eandidates has 
further inereased interest in the system. Two of the eandidates, r Ceti e and f, are near the edges 
of the theoretieal HZ. These worlds eould potentially be high priority eandidates for observations 
by a future direet deteetion mission looking for biosignatures. A more thorough analysis of the 
system is therefore warranted. We speeifieally examine the eontributions of ehemieal eomposition 
to the likelihood of e and f being habitable and having deteetable biosignatures. We examine stellar 
models with the measured eomposition and the aeeompanying habitable zone evolution. We also 
diseuss the geophysieal implieations of a high Mg/Si in planets. Although Mg/Si in r Ceti needs 
to be further eonstrained, it is elearly larger than in the Sun and Earth. 

Two main eonsiderations turn out to be important for the ease of r Ceti. Eirst, the star has a 
very high Mg/Si=l .78. The Mg/Si ratio is eloser to the boundary at whieh a planetary mantle would 
be expeeted to transition from a olivine/pyroxene mineralogy to an olivine- dominated mineralogy 
at shallow depths and a ferroperielase-rieh lower mantle. The viseosity of the mantle eould be 
mueh lower than earth. Coupled with an enhaneed abundanee of aetinide elements eontributing 
to radiogenie heating of the interior, planets around r Ceti may have mueh more vigorous and 
long-lived mantle eonveetion, whieh would in turn affeet the aetivity of surfaee proeesses and 
(bio)geoehemioal eyeling. Eurthermore, the surfaee roek eomposition may be very different from 
that of Earth as the erust may be formed by partial melting of the mantle. This eould be an essential 
point when eonsidering the habitability and deteetability of a planet 7.6 Gy old. While this eould 
benefit the habitability of t Ceti’s planets, it eould interfere with the DI of a younger planet. Easter 
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Table 2: Complete line list 


Line (A 

Element 

Lower Excitation Potential 

log (Oscillator Strength) 

EQW solar 

EQW Tau Ceti 

5052.15 

c 

7.68 

-1.3 

37.8 

8.29 

5380.32 

c 

7.68 

-1.61 

21.9 

4.65 

6156.8 

o 

10.74 

-0.43 

4.1 

1.54 

6363.79 

o 

0 

-9.72 

4.9 

1.65 

6154.23 

Na 

2.1 

-1.53 

39.8 

23.91 

6160.75 

Na 

2.1 

-1.23 

58.4 

43.52 

5711.09 

Mg 

4.35 

-1.83 

104.1 

113.7 

6318.72 

Mg 

5.11 

-1.99 

38.2 

51.6 

6696.03 

A1 

3.14 

-1.58 

38.1 

34.94 

6698.67 

A1 

3.14 

-1.95 

21.9 

19.3 

5690.43 

Si 

4.93 

-1.77 

52.6 

37.24 

5793.08 

Si 

4.93 

-2.06 

48.2 

28.3 

6125.03 

Si 

5.61 

-1.51 

34 

18.67 

6142.48 

Si 

5.62 

-1.54 

36.3 

23.49 

6145.01 

Si 

5.62 

-1.36 

41.6 

23.68 

6155.13 

Si 

5.62 

-0.78 

81.4 

55.65 

6244.48 

Si 

5.61 

-1.36 

42.4 

28.69 

6721.86 

Si 

5.86 

-0.94 

47.8 

28.65 


“Table 2 is published in its entirety in the electronic edition of ApJ, A portion is shown here for guidance regarding its 
form and content. 


Table 3: Habitable Zone Properties 


Conservative Optimistie 
AU 

Currently Habitable 

0.73 - 1.28 

0.55 - 1.32 

Continuously Habitable Range 


0.84- 1.1 

4 Gy Habitable Range 

0.67 - 1.52 

0.5- 1.6 

Habitable < -2Gy to Present 

0.73- 1.18 

0.55-1.25 
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geochemical cycling could impede the buildup of biologically produced non-equilibrium chemical 
species in the planet’s atmosphere. Time dependent analyses of planetary evolution are required. 

Second, the time and composition dependent modeling of the star provides a better estimate of 
not only the planets’ habitability potential now, but their habitability history. Given the HZ models 
used herein, planet e is very near the edge of the HZ for very optimistic assumptions. For more 
conservative cases it lies interior to the HZ at all times. Planet f highlights the time-dependent 
nature of detectability. The planet is near the edge of the predicted HZ. For the most optimistic 
case it has only recently entered, while for conservative cases it is still just outside. Even though 
it will likely see more than 7 Gy of continued habitability as the star ages, as far as we can trust 
existing climate models it has been in the HZ at most much less than 1 Gy. Earth required ~2 Gy 
before the biota modified the atmosphere sufficiently to possibly detectable from another stellar 
system. Even if f’s habitability potential is high, its detectability index is low. 

In conclusion, the known candidate planets around r Ceti are not as promising targets for find¬ 
ing life as they appear at first analysis. They do, however, show the power of full characterization 
of exoplanetary systems, in particular their chemical composition, for choosing targets for future 
missions to detect life. 
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Fig. 1.— Inner and outer radii of the HZ of r Ceti. Solid lines show the inner edge of the HZ for, 
in order of distanee from the star, the Reeent Venus, Runaway Greenhouse, and Moist Greenhouse 
eases. The outer edge is shown by the dotted lines for the Maximum Greenhouse and Early Mars 
eases. The orbits of the putative planets are indieated by the horizontal dashed lines. The vertieal 
line is the best-fit age of the system at 7.63 Gy with the age uneertainty indieated by the gray box. 









